Toward Full Spatio-Temporal Control on the Nanoscale 
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We introduce an approach to implement full coherent control on nanometer length scales. It is 
based on spatio-temporal modulation of the surface plasmon polariton (SPP) fields at the thick edge 
of a nanowedge. The SPP wavepackets propagating toward the sharp edge of this nanowedge are 
compressed and adiabatically concentrated at a nanofocus, forming an ultrashort pulse of local fields. 
The one-dimensional spatial profile and temporal waveform of this pulse are completely coherently 
controlled. 
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Two novel areas of optics have recently attracted a 
great deal of attention: nanooptics and ultrafast op- 
tics. One of the most rapidly developing directions in 
ultrafast optics is quantum control, in which coherent 
superpositions of quantum states are created by exci- 
tation radiation to control the quantum dynamics and 
outcomesi 2 -^. Of special interest are coherently con- 
trolled ultrafast phenomena on the nanoscale where the 
phase of the excitation waveform along with its polariza- 
tion provides a functional degree of freedom to control 
nanoscale distribution of energy^i 7 -^^. Spatiotempo- 
ral pulse shaping permits one to generate dynamically 
predefined waveforms modulated both in frequency and 
in space to focus ultrafast pulses in the required micro- 
scopic spatial and femtosecond temporal domain o 11 ! 12 . 

In this Letter, we propose and theoretically develop a 
method of full coherent control on the nanoscale where a 
spatiotemporally modulated pulse is launched in a graded 
nanostructured system. Its propagation and adiabatic 
concentration provide a possibility to focus the optical 
energy in nanoscale spatial and femtosecond temporal re- 
gions. The idea of adiabatic concentratio n 13 ! 14 (see also 
Ref. [HI) is based on adiabatic following by the propa- 
gating surface plasmon-polariton (SPP) wave of a plas- 
monic waveguide, where the phase and group velocities 
decrease toward a limit at which the propagating SPP 
wave is adiabatically transformed into a standing surface 
plasmon (SP) mode. This effect has been further devel- 
oped theoretically 16 ^ and observed experimentally. 18 

To illustrate the idea of this full coherent control, con- 
sider first the adiabatic concentration of a plane SPP 
wave propagating along a nanowedge of silver—, as shown 
in Fig.QJa); the theory is based on the Wentzel-Kramers- 
Brillouin (WKB) or quasiclassical approximation, also 
called the eikonal approximation in optics 2 * as suggested 
in Refs. [l3lll3 . The propagation velocity of the SPP along 
such a nanowedge is asymptotically proportional to its 
thickness. Thus when a SPP approaches the sharp edge, 
it slows down and asymptotically (in the ideal limit of 
zero thickness at the apex) stops while the local fields 
are increased and nano-concentrated. 

Now consider a family of SPP rays (WKB trajecto- 
ries) propagating from the thick side of a nanowedge, 




y(\im) 



FIG. 1: (a) Illustration of adiabatic concentration of energy 
on the wedge. The distribution of local field intensity I in the 
normal plane of propagation of SPPs (the yz plane). The in- 
tensity in relative units is color coded with the color scale bar 
shown to the right, (b) Trajectories of SPP rays propagating 
from the thick to sharp edge of the wedge. The initial coor- 
dinate is coded with color. The black curves indicate lines of 
equal phase (SPP wave fronts). 

as shown in Fig. [ljb). One possibility to launch such 
SPPs is to have nanoscale inhomogeneities (nanoparti- 
cles or nanoholes) at the thick edge of the wedge. Each 
of the optical pulses from the spatiotemporal shaper is fo- 
cused on the corresponding nanoparticle, scattering from 
it and generating an SPP wave. This scattering is neces- 
sary to impart large transverse momenta on SPPs, which 
is required for their efficient transverse focusing. These 
SPP waves propagate toward the sharp edge, adiabati- 
cally slow down which increases the field amplitudes, and 
constructively interfere as they converge at the nanofo- 
cus. 

The phases of the SPPs rays (i.e., the corresponding 
wave fronts) are defined by the spatiotemporal modulator 
in such a way that the rays converge to a reconfigurably 
chosen point at the sharp edge of the nanowedge where 
they acquire equal phases. When the SPPs propagate 
along the rays, adiabatic concentration takes place be- 
cause the SPP wavelength tends to zero proportionally 
to the thickness of the wedge. This allows one to fo- 
cus optical energy at a predefined nanofocus at the sharp 
edge. The temporal structure of the generated SPPs can 
be chosen in such a way that at the nanofocus the lo- 
cal fields form an ultrashort pulse. Using spatiotemporal 
modulation of the excitation field at the thick edge, one 
can arbitrarily move the nanofocus along the sharp edge. 
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A superposition of such fields can render arbitrary spa- 
tiotemporal modulation on the sharp edge, enabling one 
to exert full control over nanoscale fields in space and 
time. 

Turning to the theory, consider a nanofilm of metal 
in an xy plane whose thickness d in the z direction is 
adiabatically changing with the coordinate-vector p = 
{x,y) in the plane of the nanofilm. Let e m = e m (Lu) be 
the dielectric permittivity of this metal nanofilm, and £d 
be the permittivity of the embedding dielectric. Because 
of the symmetry of the system, there are odd and even 
(in the normal electric field) SPPs. ft is the odd SPP that 
is a slow-propagating, controllable mode. The dispersion 
relation for this mode defining its effective index n(p) is 



1 



taiil) ( -k d(p)y / n(p) 2 - e m j = - 



(1) 

where fco = lo/c is the radiation wave vector in vacuum. 

Let r be a unit tangential vector to the SPP trajectory 
(ray), ft obeys a conventional equation of ray optics^ 



dr 

n— = Vn — t (tvb) 
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where I is the length along the ray and V = d/dp. 

Now let us consider a nanofilm shaped as a nanowedge 
as in Fig. [Tfb). In such a case, n — n(y), and these 
trajectory equations simplify as 
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From these, it follows that n x = r x n — const. The 
SPP wave vector, related to its momentum, is k(p) = 
kon(p)T; this is the conservation of k x (the transverse 
momentum). This allows one to obtain a closed solu- 
tion for the ray. The tangent equation for the ray is 
dx/dy — t x /t v , where T y — yl — n 2 /n 2 . From this, we 
get an explicit SPP trajectory (ray) equation as 



x — Xq 
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dy' , 
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where po = (xo,ya) is the focal point where rays with 
any n x converge. To find the trajectories, as n{y) we use 
the real part of effective index (fTJ), as WKB suggests. 

When the local thickness of the wedge is subwavelength 
(k a d <C I), the form of these trajectories can be found 
analytically. Under these conditions, dispersion relation 
(TT|) has an asymptotic solution 



n = — , n a = In ■ 

k d e m + Sd 



(5) 



Substituting this into Eq. (j4]), we obtain explicit equa- 
tions of trajectories, 



x — xq 



(6) 
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FIG. 2: (a) Phase (real part of eikonal $) acquired by a SPP 
ray propagating between a point with coordinate x on the 
thick edge and the nanofocus, displayed as a function of x. 
The rays differ by frequencies that are color coded by the 
vertical bar. (b) The same as (a) but for extinction of the ray 
(Im$). 



where h a = n a /(ko tan0), and tan6> is the slope of the 
wedge. Thus, each SPP ray is a segment o f a circle whose 
center is at a point given by x = xq + y/ (n a /n x ) 2 — y$ 
and y = 0. This analytical result is in agreement with 
Fig. [T] (b). If the nanofocus is at the sharp edge, i.e., 
yo = 0, then these circles do not intersect but touch and 
are tangent to each other at the nanofocus point. 

As an example we consider a silver— nanowedge illus- 
trated in Fig. [T] (b) whose maximum thickness is d m = 30 
nm and whose length (in the y direction) is L = 5 /im. 
Trajectories calculated from Eq. ^ for fku = 2.5 eV are 
shown by lines (color used only to guide eye); the nanofo- 
cus is indicated by a bold red dot. The different trajec- 
tories correspond to different values of n x in the range 
< n x < n(L). In contrast to focusing by a conven- 
tional lens, the SPP rays are progressively bent toward 
the wedge slope direction. 

The eikonal is found as an integral along the ray 
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Consider rays emitted from the nanofocus [Fig. [T] (b)]. 
Computed from this equation for frequencies in the visi- 
ble range, the phases of the SPPs at the thick edge of the 
wedge (for y = L) are shown in Fig. [2] (a) as functions 
of the coordinate x along the thick edge. The colors of 
the rays correspond to the visual perception of the ray 
frequencies. The gained phase dramatically increases to- 
ward the blue spectral region, exhibiting a strong disper- 
sion. The extinction for most of the frequencies except 
for the blue edge, displayed in Fig. [2] (b), is not high. 

Now consider the evolution of the field intensity along 
a SPP ray. For certainty, let SPPs propagate along the 
corresponding rays from the thick edge of the wedge to- 
ward the nanofocus as shown in Fig. [T] (b). In the pro- 
cess of such propagation, there will be concentration of 
the SPP energy in all three directions (3d nanofocusing) . 
This phenomenon differs dramatically from what occurs 
in conventional photonic ray optics. 

To describe this nanofocusing, it is convenient to intro- 
duce an orthogonal system of ray coordinates whose unit 
vectors are r (along the ray), r) — {—t v ,t x ) (at the sur- 



face normal to the ray), and e z (normal to the surface). 
The concentration along the ray (in the r direction) oc- 
curs because the group velocity v g = [d(kon) / du>] of 
SPP asymptotically tends to zero (for the antisymmetric 
mode) for fco^ -* as w 5 = vo g d where vq 9 — constj 13 ! 14 
This contributes a factor Au = 1/ ' yVg (d) to the ampli- 
tude of an SPP wave. 

The compression of a SPP wave in the e z (vertical) 

direction is given by a factor of A z = I Wdz j , 

where W is the energy density of the mode. Substituting 
a standard expression 20 for W , one obtains explicitly 



A z = ^ exp (Re /c d d) | 
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where n m = kgy/n — e m and n d — ko^/n — Ed- 

To obtain the compression factor A± for the r] direc- 
tion), we consider conservation of energy along the beam 
of rays corresponding to slightly different values of n x . 
Dividing this constant energy flux by the thickness of 
this beam in the r\ direction, we arrive at 



A i = 
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The ray amplitude thus contains the total factor which 
describes the 3d adiabatic compression: A = AnA±A z . 

Now consider the problem of coherent control. The 
goal is to excite a spatiotemporal waveform at the thick 
edge of the wedge in such a way that the propagating 
SPP rays converge at an arbitrary nanofocus at the sharp 
edge where an ultrashort pulse is formed. To solve this 
problem, we use the idea of back-propagation or time- 
reversal<2i We generate rays at the nanofocus as an ultra- 
short pulse containing just several oscillations of the light 
field. Propagating these rays, we find amplitudes and 
phases of the fields at the thick edge at each frequency 
as given by the eikonal ^(p). Then we complex con- 
jugate the amplitudes of frequency components, which 
corresponds to the time reversal. We also multiply these 
amplitudes by exp(2Im $) which pre-compensates for the 
losses. This provides the required phase and amplitude 
modulation at the thick edge of the wedge. 

We show an example of such calculations in Fig. [3] 
Panel (a) displays the trajectories of SPPs calculated ac- 
cording to Eq. The trajectories for different frequen- 
cies are displayed by colors corresponding to their visual 
perception. There is a very significant spectral disper- 
sion: trajectories with higher frequencies are much more 
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FIG. 3: (a) Trajectories (rays) of SPP packets propagating 
from the thick edge to the nanofocus displayed in the xy 
plane of the wedge. The frequencies of the individual rays 
in a packet are indicated by color as coded by the bar at 
the top. (b)-(d) Spatiotemporal modulation of the excitation 
pulses at the thick edge of the wedge required for nanofocus- 
ing. The temporal dependencies (waveforms) of the electric 
field for the phase-modulated pulses for three points at the 
thick edge boundary: two extreme points and one at the cen- 
ter, as indicated, aligned with the corresponding x points at 
panel (a), (e) The three excitation pulses of panels (b)-(d) (as 
shown by their colors), superimposed to elucidate the phase 
shifts, delays, and shape changes between these pulses. The 
resulting ultrashort pulse at the nanofocus is shown by the 
black line. The scale of the electric fields is arbitrary but 
consistent throughout the figure. 



curved. The spatial-frequency modulation that we have 
found succeeds in bringing all these rays (with different 
frequencies and emitted at different x points) to the same 
nanofocus at the sharp edge. 

The required waveforms at different x points of the 
thick edge of the wedge are shown in Fig.[3](b)-(d) where 
the corresponding longitudinal electric fields are shown. 
The waves emitted at large x, i.e., at points more distant 
from the nanofocus, should be emitted significantly ear- 
lier to pre-compensate for the longer propagation times. 
They should also have different amplitudes due to the 
differences in A. Finally, there is clearly a negative chirp 
(gradual decrease of frequency with time). This is due 
to the fact that the higher frequency components prop- 
agate more slowly and therefore must be emitted earlier 
to form a coherent ultrashort pulse at the nanofocus. 

In Fig. [3] (e) we display together all three of the rep- 
resentative waveforms at the thick edge to demonstrate 
their relative amplitudes and positions in time. The pulse 
at the extreme point in x (shown by blue) has the longest 
way to propagate and therefore is the most advanced in 
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time. The pulse in the middle point (shown by green) is 
intermediate, and the pulse at the center (x = 0, shown 
by red) is last. One can notice also a counterintuitive fea- 
ture: the waves propagating over longer trajectories are 
smaller in amplitude though one may expect the oppo- 
site to compensate for the larger losses. The explanation 
of this fact is that the losses are actually insignificant for 
the frequencies present in these waveforms. What deter- 
mines the relative magnitudes of these waveforms is the 
coefficient A± of the transverse concentration, which is 
much higher for the peripheral trajectories than for the 
central ones. 

Figure [3] (e) also shows the resulting ultrashort pulse 
in the nanofocus. This is a transform-limited, Gaussian 
pulse. The propagation along the rays completely com- 
pensates the initial phase and amplitude modulation, ex- 
actly as intended. As a result, the corresponding electric 
field of the waveform is increased by a factor of 100. Tak- 
ing the other component of the electric field and the mag- 
netic field into account, the corresponding increase of the 
energy density is by a factor ~ 10 4 with respect to that 
of the SPPs at the thick edge. 

Consider the efficiency of the energy transfer to the 
nanoscale. This is primarily determined by the cross sec- 
tion <J SPP for scattering of photons into SPPs. For in- 
stance, for a metal sphere of radius R at the surface of 
the wedge, one can obtain an estimate cr SPP ~ R°/(d% n A), 
where A is the reduced photon wavelength. Setting 
R ~ d m , we estimate a SPP ~ 3 nm 2 . Assuming opti- 
cal focusing into a spot of ~ 300 nm radius, this yields 
the energy efficiency of conversion to the nanoscale of 
~ 10~ 3 . Taking into account the adiabatic concentration 
of energy by a factor of 10 4 , the optical field intensity 
at the nanofocus is enhanced by one order of magnitude 
with respect to that of the incoming optical wave. 

The criterion of applicability of the WKB approxima- 



tion is <9fc -1 /dy <C 1. Substituting k = k^n and Eq. 
([5]), we obtain a condition d m /(n a L) << 1. This condi- 
tion is satisfied everywhere including the nanofocus since 
n a ~ 1 and d m *C L for adiabatic grading. The min- 
imum possible size of the wavepacket at the nanofocus 
in the direction of propagation, Ax, is limited by the lo- 
cal SPP wavelength: Ax ~ 2ir/k w 2irdf/n a , where df 
is the wedge thickness at the nanofocus. The minimum 
transverse size a (waist) of the SPP beam at the nanofo- 
cus can be calculated as the radius of the first Fresnel 
zone: a = Tr/k x > Tr/(kon x )- Because n x is constant 
along a trajectory, one can substitute its value at the 
thick edge (the launch site), where from Eq. we ob- 
tain n x w n = n a /d m . This results in a w ird m /n a \ thus 
a is on order of the maximum thickness of the wedge, 
which is assumed also to be on the nanoscale. 

To briefly conclude, we have proposed and theoretically 
investigated an approach to full coherent control of spa- 
tiotemporal energy localization on the nanoscale. From 
the thick edge of a plasmonic metal nanowedge, SPPs 
are launched, whose phases and amplitudes are indepen- 
dently modulated for each constituent frequency of the 
spectrum and at each spatial point of the excitation. This 
pre-modulates the departing SPP wave packets in such a 
way that they reach the required point at the sharp edge 
of the nanowedge in phase, with equal amplitudes form- 
ing a nanofocus where an ultrashort pulse with required 
temporal shape is generated. This system constitutes 
a "nanoplasmonic portal" connecting the incident light 
field, whose features are shaped on the microscale, with 
the required point or features at the nanoscale. 
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